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We present first principles calculations of the atomic and electronic structure of electron-doped 
LaOFeAs. We find that whereas the undoped compound has an antiferromagnetic arrangement of 
magnetic moments at the Fe atoms, the doped system becomes non magnetic at a critical electron 
concentration. We have studied the electron-phonon interaction in the doped paramagnetic phase. 
For the Aig phonon, the separation between the As and Fe planes induces a non-collinear arrange- 
ment of the Fe magnetic moments. This arrangement is anti parallel for interactions mediated by 
As, and perpendicular for Fe-Fe direct interactions, thus avoiding frustration. This coupling of 
magnetism with vibrations induces anharmonicities and an electron-phonon interaction much larger 
than in the pure paramagnetic case. We propose that such enhanced interactions play an essential 
role in superconducting compounds close to an atiferromagnetic phase transition. 

PACS numbers: 71.18.-b, 71.38.-k, 74.20-z, 74.70.-b 



Very recently, Kamihara and co-workers [J, Q have 
reported a new family of Fe based compounds that are 
superconducting when doped. All these ROFeAs (R=La, 
Sm, etc) materials are formed by FeAs layers separated 
by insulating rare-earth-oxide layers. They are supercon- 
ducting when doped with electrons by substituting O for 
F, whereas the parent compounds are antiferromagnetic 
(for a recent review see Q)- The superconducting crit- 
ical temperature Tc is as high as 50 K f4]. The phase 
diagram has been well established experimentally both 
as a function of doping ^, and pressure [3] . The vari- 
ation and ordering of the antiferromagnetic phase with 
temperature has been measured with neutron diffraction 
Q. The phonon modes have been determined by Ra- 
man spectroscopy 9]. Photocmission experiments, and 
first principles calculations, have revealed a Fermi surface 
with two electron pockets and two hole pockets (loj . 

The precise mechanism responsible for the supercon- 
ductivity has not been yet established. Theoretical calcu- 
lations [11] rule out electron-phonon interaction to be the 
only pairing mechanism. Electron-phonon interaction in 
the paramagnetic phase can only account for a maxi- 
mum Tc of 0.8 K [lT|. Other calculations [H also rule 
out phonon mediated superconductivity. The magnetic 
structure of the parent compounds is well established ex- 
perimentally although, from the theoretical point of view, 
there are some discrepancies, between different calcula- 
tions and with the experimental results, in the magnitude 
of the magnetic moment at the Fe atoms (l3i] . 

In this work we present first principles calculations to 
analyze in more detail the electron-phonon interaction in 
the doped paramagnetic phase. The Density Functional 
[fl [isl calculations are performed using the SIESTA 
code lid [itI which uses localized orbitals as basis func- 
tions [l8|. We use non-local norm conserving pseudopo- 
tentials and a Local Density Approximation (LDA) for 



the exchange and correlation functional. The calcula- 
tions are performed with stringent criteria in the elec- 
tronic structure convergence (down to 10^^ in the density 
matrix), Brillouin zone sampling (up to 18000 fc-points), 
real space grid (energy cut-off of 500 Ryd) and equilib- 
rium geometry (residual forces lower than 10^^ eV/A). 
Due to the rapid variation of the density of states at the 
Fermi level, we used a special smearing method [l9| . 

We first calculate the variation of the total energy with 
the lattice parameter for various possible magnetic con- 
figurations. In Figure [1] we show the results for two pos- 
sible antiferromagnetic configurations in the Fe atoms: 
direct antiferromagnetic ordering between nearest neigh- 
bor Fe atoms (AFO) and antiferromagnetic Fe-Fe inter- 
action mediated by As (AFl). We find that, at the cal- 
culated equilibrium lattice constant (3.97A), the ground 
state has AFl magnetic order. Its energy difference with 
AFO (which is non magnetic at that lattice constant) is 
26 meV per formula unit. Ferromagnetism and other 
non-collinear magnetic configurations are found to be un- 
stable. The AFl magnetic moments are 1.3 fiB per Fe 
atom, in reasonable agreement with previous all-electron 
density functional calculations 2^, 21 1. This calculated 
value is larger than the 0.36 /is/atom found in neutron 
diffraction experiments Q , possibly due to neglected spin 
fiuctuations 
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In order to simulate the doped phase we have added 
electrons to the system while keeping the charge neu- 
trality by adding a uniform positive charge background. 
As shown in Figure [21 the magnetic moments decrease 
with electron doping, and they disappear above a criti- 
cal concentration. This critical value (0.26 electrons per 
formula unit) is larger than the experimental one, prob- 
ably due to the overestimation of the magnetic moment 
in the undoped phase. From the structural point of view, 
the main effect of electron doping is to decrease the dis- 
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FIG. 1: (Color online) Variation of total energy (left scale) 
and magnetic moment of Fe atoms (right scale) with the lat- 
tice constant a in two different antiferromagnetic arrange- 
ments of LaOFeAs: (a) direct Fe-Fe antiferromagnetic order 
(AFO); and (b) As- mediated antiferromagnetic order (AFl). 
The zero of energy is at the AFl minimum. The orienta- 
tions of the Fe magnetic moments in each configuration are 
indicated in the insets, where squares and triangles stand for 
As atoms above and below the plane of Fe atoms. The total 
energy is for a (\/2a x \/2a x c) supercell. 



tance between the As and Fe planes. For instance, in 
the paramagnetic phase, with 0.275 extra electrons, the 
compression is 0.02 A. This compression, added to the 
filling of the Fe d levels and the shift of the Fermi energy 
to a region with lower density of states, is responsible for 
the disappearance of magnetism. Figure [2jb) shows that 
the extra c harg e goes mainly to the As-Fe-As layer, as 
anticipated |22l |. 

To address the electron-phonon interaction in this sys- 
tem, we have considered the symmetric out-of-plane FeAs 
Aig mode at fcn = 0, in which the Fe atoms remain fixed 
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whereas the As atoms move perpendicularly to the FeAs 
layers, expanding and compressing the Fe-As bonds. Fig- 
ure [3] shows the calculated total energy and magnetic 
moment, as a function of the expansion (positive u) or 
compression (negative u) of the Fe-As interplanar dis- 
tance, for an extra concentration of 0.275 electrons per 
formula unit. Several points are worth mentioning: i) 
Expansion of the Fe-As distance induces magnetism in 
the Fe atoms, ii) In this case the arrangement of the mag- 
netic moments at the Fe atoms is non-collinear, such that 
the As-mediated interaction remains antiferromagnetic, 
while there is no longer magnetic frustration for the direct 
Fe-Fe interaction, as in the AFl order, in) Compression 





1 .0 - 




1 4 - 




1.2- 






c 




a 




o 


0.8 - 


E 




o 


0.6 - 


% 




cn 


0.4 - 


to 






0.2 - 




0.0- 




0.12- 




0.10- 



S 0.08 H 
S 0.06 



a> 0.04 

LU 

0.02- 



0.00 



(b) 



^1- 





— ■— Fe 






—•—As 






-*-La 






-T-0 






• 






■ 


,A 








A' 






T' 







0.00 0.04 0.08 0.12 0.16 0.20 0.24 0.28 0.32 
Electron concentration per ctiemical unit cell 

FIG. 2: (Color online) Calculated effects of electron doping in 
the AFl magnetic configuration of LaOFeAs (see text), (a) 
Magnetic moment per Fe atom, (b) Distribution of the extra 
charge among the different atoms. 



of the Fe-As planes does not induce magnetism, and this 
asymmetry makes the vibration strongly anharmonic. 
The energy per formula unit can be fitted by the fourth 
order polynomial E{u) = fiwix^ — 0.107a;'^ — 0.064x''), 
where hw = 29.0 meV is the calculated phonon energy, 
X = u/uq and uq — 0.062 A. Thus, the restoring force 
is -1-0.91 and -0.68 eV/A for u =-0.06 and -f0.06 A, re- 
spectively, iv) The magnetic moments appear even for 
small phonon amplitudes, such that E{u) < huj {u < uq). 
v) For positive u, the non coUinear solution is less than 
~ 1 meV lower than the AFl configuration, which has 
very similar magnetic moments, vi) A similar, although 
weaker, effect has been obtained in the calculation of the 
asymmetric A2u mode. 

Next we study how the Aig vibration affects the elec- 
tronic structure in detail. Figure H] shows the calculated 
band structure, with an electron doping of 0.275, for dif- 
ferent magnetic states and vibration displacements. The 
electron and hole pockets appear both at T, since the 
calculations were performed with a {\/2a x V2a x c) su- 
percell, with the primitive M point folded to F. The 
band structure is not much perturbed in the paramag- 
netic state (Figure |4fb)): the maximum deformation po- 
tential is ~ 2 eV/A and it occurs for states far from the 
Fermi level Ep. Around Ep this value is much smaller, 
in agreement with Boeri et al Uj. In contrast, the elec- 
tronic structure is strongly perturbed in both the AFl 
and non-collinear antiferromagnetic states, which remove 
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FIG. 3: (Color online) Variation of energy and magnetic mo- 
ment with the amplitude u of the Aij-phonon at F (vibration 
of As atoms perpendicularly to the Fe plane) for a doping 
of 0.275 electrons per formula unit, (a) Energy per formula 
unit in the paramagnetic and non-collinear magnetic configu- 
rations, (b) Non-collinear magnetic moment of the Fe atoms. 
Their orientations are indicated in the inset, where squares 
and triangles stand for As atoms above and below the Fe 
plane. 



several band degeneracies (Figures IDJc) and (d)). In addi- 
tion, antiferromagnetism removes states from about 0.3 
eV below Ep, while a new peak in the density of states 
appears close to Ep (Figure O. This new peak is due to 
to the flattening of the electron pocket bands at the F 
point, right at Ep- The peak size increases with u, and 
it crosses Ep for large u (see Figure [5]). The analysis of 
the partial density of states reveals that the states asso- 
ciated with this peak are fully spin polarized and entirely 
localized in the Fe atoms. 

Finally, we suggest a possible general mechanism for 
BCS-like high-temperature superconductivity in antifer- 
romagnetic compounds. It is based on driving them 
very close to a critical point or line of the antiferro-to- 
paramagnetic transition, in the pressure-doping phase di- 
agram (we assume here T = 0). Under these conditions, 
the lattice distorsion of some phonons will switch anti- 
ferromagnetism on and off, thus enhancing dramatically 
their deformation potential and the A factor [23| . In prac- 
tice, chemical doping is inhomogcneous in the nanoscale, 
what implies that the superconducting band gap will vary 
spatially, or even that patches of normal and supercon- 
ducting regions will coexist [24|] . 

In summary, we have shown that the electron-phonon 
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FIG. 4: Electronic band structure for LaOFeAs, doped with 
0.275 electrons per formula unit, for different magnetic con- 
figurations (Paramagnetic, AFl, and AF non collinear, see 
text) and vibration displacements u of the A\g phonon at P. 
The Brillouin zone is for the {\/2a x %/2a x c) supercell. The 
zero of energy is at the Fermi level. 
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FIG. 5: (Color online) Densities of states of the stable non- 
coUinear antiferromagnetic sate, for different vibration dis- 
placements u of the ^ig-phonon at F. For clarity, the curves 
have been shifted and smoothed with a gaussian broadening 
of 0.1 eV. Lower to upper curves correspond to increasing 
values of it. 



interaction in paramagnetic electron-doped LaOFeAs is 
more complex and, probably, dramatically larger than 
previously predicted. The complexity arises from the 
vicinity of the system to magnetism and, in particular, 
with the abrupt appearance of phonon-induced magnetic 
moments at the Fe atoms, in an antiferromagnetic As- 
mediated configuration. The asymmetry of magnetism, 
with respect to the compression or expansion of the Fe- 
As bonds, makes the vibrations anharmonic. These re- 
sults suggest that electron-phonon coupling must be care- 
fully revised, before ruling out its connection with high- 
temperature superconductivity, in layered compounds 
close to an antiferromagnetic transition. Calculations of 
the complex electron-phonon interaction are under way, 
in order to address quantitatively this proposed mecha- 
nism. 
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